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h ' ABSTRACT 

> ; 

■<sj- , We propose to investigate flavour changing neutral currents in the c — > wy transition 

through the measurement of the difference between T(D° — ► p°7) and T(D° — > uj). This is 
based on the observation that D° — > (dd)j is due to long distance physics while D° — > (uujj 
arises from the c — > wy transition. The effect of p — uj mixing is included. A difference in 
the decay widths of more than about 30% would be indicative of new physics. 
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Despite the remarkable success of the standard model (SM), it is generally believed that 
. £h ! this is an effective theory at present energies. The lack of explanation for many of the salient 
features of SM suggests that one must look for its extension. At present, there is no clear 
picture of the form of the SM extension and the search for physics beyond the SM proceeds 
along many alternative paths. 

As the couplings of Z°, photon and Higgs boson are flavour diagonal, flavour changing 
neutral current (FCNC) transitions are rare in SM since they can arise at loop level only. 
Moreover, FCNC transitions in the up-quark sector (like c — > wy ) are much rarer than 
those in the down-quark sector (like b —>■ 57) as a result of the CKM matrix elements and 
masses involved ]I] . Accordingly, they could provide an appropriate ground for the search of 
new physics effects 0. In fact, rare charm decays have been frequently marked as possible 
sources for the discovery of new physics 0, |3j, in view of the smallness of the short distance 
SM contributions of these processes. 

Hadronic decays like D — > V'y [1, 4-9] and D — > Vl + l~ [9-12] have been considered with 
the aim of investigating the size of the short distance (SD) c — > wy, c — > ul + l~ contribution 
to the decay amplitude. However, it turns out that these decays are largely dominated by 
long distance (LD) contributions, rendering them inappropriate for detecting deviations from 
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the SM values of the basic c — > wy(l + l ) transitions. An exception is B c — > B*^ since in this 



decay both the SD and LD contributions to the branching ratio are in the 10 range |L3 . 

In the present letter we suggest a new possibility for the search of the c — > wy transition, 
from the determination of the difference in size of the partial decay widths D° — > p°7 and 
D° — > uj. Our method is particularly suitable for detecting an enhancement of the c — ► wy 
amplitude coming from new physics if it increases this amplitude by at least a factor 3-4. We 
remark that enhancements of up to a factor 100 have been noted in certain non-minimalQ 
realizations of supersymmetry [0 . 

Next we note that in the radiative decays of D° mesons, the (dd) r j final state arises 
primarily as a result of nonleptonic W- exchange cu — > dd, being therefore an outcome 
of LD physics (Fig la), while the (uu)^ final state is mainly due to the electro-magnetic 
penguin c — > wy transition (Fig. lb). The {uu) r ) final state receives also a small contribution 
form the long distance penguin mechanism illustrated in Fig. lc. The decays D° — > (p, oj) r y 
have been treated in great detail in Refs. [0, ||, |9| and it is shown that the LD contributions 
completely overshadow the SD contribution in SM. Although, there is a strong enhancement 
at the two - loop level |TJj] of the c — > wy amplitude |TJ] as a result of gluonic corrections, 
the QCD - corrected SM amplitude |L4] is still only a few percent of the LD one ||. 



As usual, we define the isospin eigenstates 

u (i=o) = J_( Uu + ^d) , p( /=1 ) = - dd) . (1) 

v 2 V 2 

Therefore, in the absence of the uj — p mixing and the LD penguin contribution, the LD rates 
for D° — > p7 and D° — > wy are equal in the SM; a difference must be due to SD effects. This 
is the basis for our approach. 

The main decay modes of the physical p, uo states are 2n and 3ir respectively. It is well 
known that the cu^ I=0 ^ and p( /=1 ) states can mix and u is known to decay also to 2ti states, 
with a branching ratio of Br(u — > 7r + 7r~) = (2.21 ±0.30)% ]nj. The physical states denoted 



as p and uj are related to the states defined in Eq. ([J) by ]T7| 



p = p (/=D _ e o,(/=o) } u = w (/=o) + e p (/=i) j (2 ) 

where 



n 2 n 2 



ml - m 2 m 2 - m 2 + im p T p - im w T u 



(3) 



with complex masses rh = m + imT/2 and with real m and T. This leads to the uj — > 27r 
amplitude 

^(w ^ 2tt) = — j-Si- —— A(p - 2tt), (4) 

from which the mixing parameter is determined to be 

n 2 = -(4.0 ± 0.4) x 10- 3 GeV 2 (5) 



^^Non- minimality here means that the universality of the soft breaking terms is not imposed or that 
additional super-fields are added to MSSM. 
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where the values of m p , and T p are taken from [16]. The minus sign in (H) is obtained 
from a detailed analysis of the e + e _ — > 7r + 7r~ amplitude [[17]]. The corresponding value of e 

(S)is 

e = -0.0061 +i 0.036, (6) 
and the physical states can be written as 

p= -j=(l-e)uu+ -j=(-l - e)dd u = -j={1 + e)uu + -j=(l - e)dd . (7) 

The general form of a D — > V^7 transition amplitude is written as usual in terms of parity 
conserving Ape and parity violating Apy amplitudes [7|, |8| 

+a PV [(e^ ■ g)(e; • p v ) - (q ■ Pv )(e v ■ e*)]}. 
We decompose the Ape and Apy amplitudes according to their final state content 

A PC = A PC {uu) + Ap C {dd), A PV = A PV (uu) + A PV (dd) . (9) 
At this point we define the new quantities rjpc, Vpv 

A PC (uu) A PV (uu) 
VPC ~ A^ddj ' VPV ~ AMddj ' (10) 

which, as we explain later on, are of the order of a few percent in the standard model, and 
the amplitudes are rewritten as 

Apc/pv{D° -> p° 7 ) = -^=A PC /Pv(dd) [(1 - e) VPC/PV + (-1 - e)] , 
Apc/pv(D -> loj) = -^=A PC/PV (dd) [(1 + e)7]p C/ p V + (1 - e)]. (11) 
The decay width is given by 

r(D - v-r) = ^^^n\ A pc\ 2 + \Apv\ 2 ) • (12) 

In order to extract the (uu)^ final state, we propose a quantity V w ~ p , defined as 



T[D° -> ur(]/(m 2 D - ml) 3 - T[D° -> p° 7 ]/(m 



2 _ m 2^3 



r[D°^o;7]/K-m2)3 ' 1 J 

The standard model values of t]pc/pv (given in (^) bellow) and e (|5|) are small and D u '~ p 
can be expanded to the first order 



|-4pc(^)| 2 (Re r/ PC - Re e) + 


ylpy((i(i)| 2 (Re rjpv - Re e) 




>lp C (rf(J)| 2 + |A PV -(rfJ)| 2 
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which leads to the near equality of the D° — > p7 and D° — > wy rates. On the other hand, 
physics beyond the standard model, which affects significantly the size of c — > wy transition, 
will lead to a splitting of the T[D° — > p°7] ~ T[D° — > CJ7] near equality. Sizeable coefficients 
can arise in such scenarios and the expansion (0) is not valid. We remark that as a 
consequence of the mixing term e 7^ (Q) there is a difference of a few percent between the 
widths of the two decays, even when r]pc and r]py are neglected. An additional difference 
arises at the hadronic level as a consequence of SU(3) flavour breaking. For example, as one 
sees from the explicit expressions of 0, the difference in couplings and masses of p and uj 
induces a ~ 10% difference between the D° — > p°y and D° — > uy widths. 

We estimate the long distance contribution A(dd) from the calculation of the decays D° — > 
p°7, wy derived by the use of the HQET and chiral Lagrangian together with factorization 
in Refs. [[F], [8]. The relative signs of different contributions to the parity conserving and 
parity violating amplitudes were left undetermined in Refs. |7], || and are taken from the 
quark models, as in ||, giving 

Apc(dd) ~ -3.7 x 1(T 9 GeV" 1 , A PV (dd) ~ 9.5 x 10~ 9 GeV" 1 . (15) 

Inserting the numerical values of A(dd) (O) and e (H) to the expression for X>" J ~ P (IB) we get 



V 



w-p , 



0.025 + 0.54 Re7]p C + 3.6 Rer] PV - 0.003 


Vpc 


2 - 0.02 


Vpv 


2 


1 + 0.27 Rerjpc + 1.8 Rer/ Py - 0.02 Im^c - 0.13 lmrj PV 


+ 0.14 


Vpc 


2 


+ 0.9 


?7pv 


2 



(16) 



where all orders in r] PC /pv are retained. This expression displays the sensitivity of the 
quantity defined by (0), on the magnitude of the short distance rate c — > wy contained 

in rjpc/pv C|10D - W e n °t e that even if we neglect r/pc/pv, there is a difference of about 2.5% 
in the decay widths due to the p — uj mixing (^|, ||), and another possible difference coming 
from SU(3) breaking; together these amount to at most 15%. 

The standard model prediction for the short distance contribution A SD (uu) is extracted 
from the calculation of the c — > wy rate at the two loop level |TJ| giving Br(c — > wy) ~ 
3 x 10~ 8 for D° decays. The matrix element (p°,uj\ua l j >1/ (l +75)c|-D°) is calculated using the 
procedure described in detail in Ref. 0, giving 

Afg(im) ~ (-1.4 x 10~ 10 - i 4.0 x 10~ 10 ) [1 ± 0.2] GeV" 1 , 

A s P %(uu) ~ (-2.3 x 10~ 10 - i 6.6 x 10~ 10 ) [1 ± 0.2] GeV" 1 . (17) 

The mechanism, shown in Fig . lc, gives rise to a long distance part of A(wu) and has been 
estimated in |7], |{| 

Aj£(uu) ~ (2.2 ± 2.2) x 10- 10 GeV" 1 , Ap^(uu) ~ (3.7 + 3.7) x 10~ 10 GeV" 1 . (18) 
This contribution is small as it is proportional to the SU(3) flavour breaking parameter 

2m 2 p 6m~ 3m^ 
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with g v defined as (V(q, e)\jy\0) = gy(q 2 )t* p . The mean value and error in ( |19D are deter- 
mined from the experimental data on V° — > e + e~ decays |16| by assuming 5v(0) ~ gvi^v) 



Given that A{uu) = A SD {uu) + A LD {uu) (|I7|, |18|), the standard model prediction for r] PC /pv 
is 

Re r/ PC - -0.02 ± 0.06 , Im r} PC ~ 0.11 ± 0.02 (20) 
Re 7]p V ~ 0.01 ± 0.04 , Im r] PV ~ -0.07 ± 0.01 

leading to the approximate equality 

Br[D° -> pj] ~ Br[D° -> cj 7 ] ~ 1.2 x 10" 6 . (21) 



The standard model prediction for the relative difference of the rates T> u p ( 13| , |16|) is 



Vffi ~ 6 ± 15 % , (22) 

where the error is dominated by the uncertainty of the ^[/(S) flavour breaking parameter, 
including (|T9"D. The strong rescattering, which can transform a dd pair to a mm pair does not 
change this result. This is because we consider decays into strong interaction eigenstates p 
and uj, which evolve as exp[— i(E — i^T)t]. 

For completeness, we consider also the case, when the p and uj states are experimentally 
identified only by the 2n and 3ti decay modes, respectively, and not by their masses and 
widths. If the two pion final states, which arise from the decay of uj, are not disentangled 
in the experiment, then one is dealing with the final states of good isospin. A final state, 
which starts out as p^ I=l \ will obtain an admixture of uj( i= ^ due to the isospin mixing (|2|), 
and vice versa. The production of 27r7 and 37r7 final states depends on time. In order to 
extract the short distance contribution, we propose to look at the number of 2n r j and 37r7, 
integrated over time 



N[D°^3w 7] _ JV[£>°->27r 7] 

7,371-2* _ {m 2 D -mlfBr(u}->Z-K) {m 2 D -m? p f 
U - N[D»->3-k 7] ' 

where 2n and 3-7T have invariant masses covered by p and uj resonances. At the first order in 
Vpc,pv and e we get in the standard model 



\A PC (dd)\ 2 Re r] PC + 


A PV (dd)\ 2 Re r] PV 




A PC {dd)\ 2 + 


A PV (dd)\ 2 



V sm * 4 '~ \r)A\ 9 — * 4± 15 % , (24) 



where we have taken into account that e{E UJ — E p )/T p <C t)pc,pv ■ 

The quantities T> w ~ p (li~3| ) and T> 3n ~ 27T (|2~3"D are particularly sensitive to new physics scenar- 
ios, which could enhance the c — > wy rate. Non-minimal realizations of the supersymmetric 
models, discussed in [E|, can enhance these quantities up to 



Similar reasoning may be useful in extracting the FCNC transition c — * ul + l~ from the 
difference in the decay rates T[D° — >p°l + l~] and T[D° — >ujI + 1~], where I denotes an electron 
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or muon. Long distance contributions to these have been studied in [RL |12j and arise via the 
mechanisms illustrated in Figs, la and lc, where the photon is replaced with the virtual 
photon decaying to leptons. The corresponding differential branching ratio has resonant 
shape in terms of the di-lepton mass m« with maximums at m\i = m p , m w and m^. The 
LD contributions, which give rise to the final state (dd)l + l~, largely cancel in the difference 
F[D° -> - T[D° -> pl + l~], which is proportional to A(uu) = A SD (uu) + A 



LD, 



uu) 



In this case, however, the SU(3) flavour cancellation in A LD (uu) is not so effective as the 
maximums at mii — m P)UJ and mu = m$ are well separated^]. As a consequence, \A LD [uu)\ is 
more than one order of magnitude larger than | A SD (uu) | and overshadows the interesting 
FCNC transition c -> ul + l~ in the difference T[D° -> -V[D° -> pl + l~]. 

To summarize, we have proposed here yet another test for physics beyond the standard 
model in the charm sector. Our test exploits the extreme smallness of the c — > wy transition 
in standard model JI3| and the near equality of the D — > p°7, D° — > wy amplitudes, obtained 
from the calculation 0, |], || of long distance contributions to these decays. This equality 
would be spoiled, if, as encountered in certain supersymmetric models [0, the c — > wy 
amplitude is enhanced by a seizable factor. Standard model calculations, including SU(3) 
breaking in the long distance contribution to the amplitudes, show that D — > p°7, wy widths 
do not differ by more than 15%. Thus, we are led to claim, conservatively, that a difference 
larger than 30% would be a "smoking gun" indication of new physics. The formalism leading 
to this conclusion has been exposed in this paper. In view of the expected ~ 10~ 6 branching 
ratio for these decays, we look forward to the experimental tests, keeping in mind that the 
present upper limits are around 10~ 4 |1 



The standard model predictions presented in [| g give Br LD {uu) - 1CT 7 and Br SD ( uu) ~ 10 . 
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(a) Long distance contribu- 
tion leading to the (dd)~f fi- 
nal state. The photon can 
be emitted from any of the 
quark lines. 



U 



u 



(b) Short distance contribu- 
tion arising from the FCN 
decay c — > wy. The hexagon 
denotes the loop diagrams 
responsible for the c — > wy 
transition. 




u 



u 



(c) Long distance contribu- 
tion leading to the (uu)^ 
final state. The box de- 
notes the action of the 
weak nonleptonic effective 
Lagrangian. 



Figure 1: Illustration of different contributios to the decays D° —> p°7 and D° — > to-y. 
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